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ABSTRACT 

The rate of terrestrial irradiation events by galactic gamma-ray bursts (GRBs) is estimated 
using recent standard-energy results. We assume that GRBs accelerate high-energy cosmic rays, 
and present results of three-dimensional simulations of cosmic rays moving in the Galactic mag- 
netic field and diffusing through pitch-angle scattering. An on-axis GRB extinction event begins 
with a powerful prompt 7-ray and neutron pulse, followed by a longer-lived phase from cosmic-ray 
protons and neutron-decay protons that diffuse towards Earth. Our results force a reinterpre- 
tation of reported ~ 10 18 eV cosmic-ray anisotropies and offer a rigorous test of the model 
where high-energy cosmic rays originate from GRBs, which will soon be tested with the Auger 
Observatory. 

Subject headings: astrobiology — cosmic rays — gamma-rays: bursts 

1. Introduction 

Observations link GRBs — brief flashes of 7-ray light emitted by sources at cosmological distances — with 
star-forming galaxies and high mass stars (van Paradijs et al. 2000). These results indicate that GRBs 
are produced by a rare type of supernova where the evolved core of a massive star collapses to a black 
hole. Normal supernovae, where the stellar core collapses to form a neutron star, are thought to produce 
nonrelativistic supernova remnant shocks that accelerate cosmic rays with energies < f 14 eV. The origin of 
higher energy cosmic rays is controversial. One possibility is that high-energy (> 10 14 eV) cosmic rays are 
accelerated by the relativistic shocks formed by GRB explosions in the Galaxy and throughout the universe 
(Vietri 1995; Waxman 1995; Dermer 2002; Wick et al. 2004), dating back to earlier suggestions of a Galactic 
origin of high-energy cosmic rays (Kulikov et al. 1969). 

If GRBs accelerate high-energy cosmic rays, then Galactic GRBs could be detected as cosmic-ray sources 
from neutron /3-decay emissions (Ioka et al. 2004). GRBs with their jets oriented toward the Earth have 
potentially lethal consequences, and may have contributed to past extinction episodes (Dar et al. 1998; 
Melott et al. 2004). The sources of the « 10 18 eV cosmic-ray excesses reported from measurements made 
with the SUGAR (Sydney University Giant Air Shower Recorder) and AGASA (Akeno Giant Air Shower 
Array) detectors (Bellido et al. (2001); Hayashida et al. (1999); see Nagano & Watson (2000) for review) 
have been proposed to result from past GRBs in the Galaxy (Biermann et al. 2004). 

To investigate these ideas, the rate of GRB events at different fluence levels is estimated, based on recent 
findings about beaming in GRBs (Frail et al. 2001; Bloom et al. 2003). A 3D propagation model is used 
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to simulate the sequence of irradiation events that occurs when a GRB jet is pointed towards Earth. The 
results suggest that a GRB jet could have produced radiations that contributed to the Ordovician extinction 
event (Melott et al. 2004). Our results are contrary to the claim that the SUGAR excess could be produced 
by a GRB in the Galaxy. If the Auger Observatory confirms the SUGAR cosmic-ray point source, then a 
model of high-energy cosmic rays from galactic GRBs is incomplete. 



2. Rate of GRB Events at Different Fluence Levels 

Most of the electromagnetic radiation from a GRB is emitted during the prompt and early afterglow 
phases on timescales of minutes to hours. We define the bolometric photon fluence <p — Sip Q with reference 
to the Solar energy fluence ip® = 1.4 x 10 6 <S ergs cm~ 2 received at Earth in one second. Significant effects 
on atmospheric chemistry through formation of nitrous oxide compounds and depletion of the ozone layer is 
found when S > 10 2 - 10 3 (Ruderman 1974; Thorsett 1995; Gehrels et al. 2003; Thomas et al. 2005), taking 
into account the very hard incident radiation spectrum of GRBs. Reprocessing of incident GRB radiation 
into biologically effective 200 - 320 nm UV radiation (i.e., with 1/e lethality) on eukaryotes occurs when 
S > 10 10 2 (Scalo & Wheeler 2002; Smith, Scalo, & Wheeler 2004). 

Achromatic beaming breaks in GRB optical/IR afterglow light curves (Stanek et al. 1999), if due to 
jetted GRBs, imply typical GRB jet opening half-angles (0j) = 0.1. Analyses (Frail et al. 2001; Bloom et 
al. 2003) show that long-duration GRBs have a standard total energy £ = 6^E lt i so /2 = 10 51 £5i ergs, with 
£51 = 1.33, a variance by a factor of 2.2, and a low E^j so population. Here E 7tiso is the apparent isotropic 
7-ray energy release inferred directly from observations. If the cone of emission from a GRB at distance 
R intercepts the line-of-sight to Earth, then the radiant fluence is given by ip = E 7 ^ so /4:TrR 2 . Thus the 
maximum sampling distance R s of a GRB with apparent isotropic 7-ray energy release E ljiso to be detected 
at the fluence level <p > ip t h — Sip Q is 



E^.jso ^ 1.1 kpc / £51 

V A ^ th ~ to/o-i) YOS/io 3 ) • [i) 

When S < 10 5 , the sampling distance to a typical GRB exceeds the « 100 pc disk scale height of 
molecular clouds and OB associations, and we can approximate the distribution of GRBs in the Galaxy by 
a uniform disk of radius Rmw — 15i?i5 kpc. The fluence size distribution of GRBs in this approximation is 
given by 

N(>S)^N GRB (-^fV(e ] ), (2) 
Kmw 

where Ngrb is the rate of GRBs in the Milky Way, and V(0j) — # 2 /2 is the probability that the Earth lies 
within the emission cone of a 2-sided jet. If GRBs follow the star formation rate history of the universe, 
then one GRB occurs every 10 5 t 5 years in the Milky Way (Wick et al. 2004; Dermer 2002), with t 5 ~ 0.1 - 
1. Thus 

Estimates (1) and (3) show that a GRB at a distance w 1 kpc with S 3> 10 2 takes place about once every 
Gyr, and more frequently if t§ = 0.1. 
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3. Cosmic Ray Propagation in the Galaxy 

We have developed a numerical model where cosmic rays move in response to a large-scale magnetic field 
that traces the spiral arm structure of the Galaxy, and diffuse through pitch-angle scattering due to magnetic 
turbulence. The magnetic field B of the Galaxy is modeled as a bisymmcrtric spiral for the Galaxy's disk, 
and a dipole magnetic field for the Galaxy's halo (Alvarez-Muhiz et al. 2002). The evolution of the particle 
momentum p = mjf3c is found by solving the Lorentz force equation dp/dt = q(i x B, where q and m are 
the particle's charge and mass, respectively, /3c is its velocity, and 7 = (1 — /3 2 ) -1 / 2 . 

Magnetic turbulence causes a particle to change its pitch angle by w tt/2 when travelling the mean- free- 
path A. The energy dependence of A is obtained by extrapolating the expression for A in a diffusion-model 
fit (Wick et al. 2004) to the measured ionic flux near the knee of the cosmic-ray spectrum (Kampert et al. 
2001) to high energies. Our approach assumes isotropic turbulence that is uniform in the disk and halo of the 
Galaxy; anisotropic turbulence is more realistic (Goldreich & Sridhar 1997), but increases the number of free 
parameters in the model. The particle's azimuth and cosine angle are randomly chosen between and 2tt 
and between \i m i n and 1, respectively, after every A[(l — n m in) I 2] 2 / cAt steps, where At is the time interval 
of each step in the numerical integration, which is set equal to a small fraction of the local gyroperiod. 

The propagation calculations are performed in the test-particle approximation, assuming that the cosmic 
rays do not affect their surrounding environment. The energy density of the cosmic-ray shell can, however, 
greatly exceed the magnetic field energy density of the ISM, in which case the protons will sweep up material, 
causing strong adiabatic losses and rapid deceleration withn a Sedov length scale of ~ pc. This forms a shock 
and a second phase of cosmic ray acceleration, which will require a hydrodynamic simulation to treat. The 
cosmic-ray neutrons decay over the entire jet volume. Comparing the ISM magnetic-field energy density 
with the energy density of the neutron-decay protons over their decay volume shows that the propagation 
of neutrons with 

> 2 x 10 7 

1 '~ (B/3^G)°- 62 (l-cos^) - 31 ' [ ' 

Thus our subsequent discussion of ~ 10 18 eV cosmic-ray neutrons is accurately described in the test particle 
approximation. 

The propagation of cosmic-ray protons and neutrons, with and without the effects of diffusive scattering, 
is illustrated in Fig. 1. Cosmic-ray neutrons travel ~ 10(i?/10 18 eV) kpc before decaying. Cosmic-ray protons 
with energies > 3 x 10 18 eV and cosmic-ray neutrons with energies > 10 18 eV escape almost directly into 
intergalactic space. Protons with energies < 3 x 10 17 eV diffusively escape from the Galaxy through a 
combination of Larmor motions and pitch-angle scatterings. 

Fig. 2 displays the cosmic-ray halo that surrounds a GRB source ~ 4 x 10 11 s after the event. The 
GRB is modeled by radially-oriented jets with 9j = 0.1. A conical shell with effective angular extent larger 
than 0.1 radians forms as a result of directly accelerated protons and neutron-decay protons with E > 10 18 
eV. The turbulent wave spectrum that is resonant with these high energy cosmic rays is poorly known. In a 
simulation with pitch-angle scattering absent, energy-dependent features and wall-like structures are formed, 
as shown in the inset to Fig. 2. 
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Fig. 1. — Trajectories of cosmic rays ejected outward from a cosmic-ray source located at 3 kpc from the 
center of the Galaxy. The spiral curves show the peak magnetic field values in the disk of the galaxy for the 
magnetic field model of Alvarez-Muhiz et al. (2002). The orange and turquoise curves show trajectories taken 
by cosmic-ray protons with Lorentz factors 7 = 10 9 and 3 x 10 9 , respectively, with diffusive scattering omitted, 
and the green and blue curves illustrate the effects on protons with respective energies when diffusive pitch- 
angle scattering is included. The red curve shows the path taken by a cosmic-ray neutron with 7 = 3 x 10 s 
and mean decay length of 3 kpc, including scattering effects on the path of the neutron-decay proton. 
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Fig. 2. — (top) Cosmic-ray halo formed 12,000 years after a GRB that took place at 3 kpc from the center of 
the Galaxy. For clarity, equal numbers of cosmic rays are injected per decade, with cosmic rays color-coded 
by energy. Cosmic-rays protons with Lorentz factors 7 in the ranges 10 7 - 10 8 , 10 s - 10 9 , 10 9 - 10 10 , and 10 10 
- 10 11 are red, green, turquoise, and purple, respectively. Cosmic-ray neutrons and neutron-decay protons 
are dark red, yellow, magenta, and dark blue in the respective energy ranges. (Bottom) Side view of the 
cosmic-ray halo when diffusive scattering is turned off (grid lines are 2 kpc in separation) . 
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4. Discussion and Summary 

About once every several hundred million years, the Earth is illuminated by prompt GRB photon 
and neutral radiations of sufficient intensity to have significant effects on the biota through erythema from 
cascade UV flux (Scalo & Wheeler 2002; Smith, Scalo, & Wheeler 2004) and destruction of the ozone layer 
catalyzed by the formation of nitric oxide and NOj, "odd nitrogen" compounds (Ruderman 1974; Thorsett 
1995; Thomas et al. 2005). The destruction of the ozone layer causes the Solar UV fluence to greatly exceed 
the GRB cascade UV fluence when the net effects of the GRB radiation on the atmosphere are considered 
(Thomas et al. 2005). Such an event might have been responsible for trilobite extinction in the Ordovician 
era through the destruction of plankton (Melott et al. 2004). 

If high-energy cosmic-ray production accompanies GRBs, then Galactic events could also have affected 
biological evolution due to DNA radiation damage by the elevated ground-level muon fluxes induced by 
w 10 17 eV (78 = 7/IO 8 = 1) cosmic-ray neutrons (Dar et al. 1998). The prompt energy fluence of neutrons 
is 

2 2 dN n exp(-r/r„) _ 4 x 10 8 T] n £ 52 cxp(-r kpc /j 8 ) -2 ,^ 

^' r ) =m ^^T4^ = T< ergs cm , (5) 

where an on-axis GRB at a distance of rk pc kpc accelerates 10 52 r) n £52 ergs of nonthermal neutrons above 
w 100 TeV into a two-sided jet with beaming factor /_2%- The prompt muon number flux due to this type 
of event is at the level 

, p ^ 4xl0 n scc6> 77„£ 52 11s ia\ 

(Gaisser 1990) for GRBs within tt/3 of the zenith. This value is above the level <p„(> 3 GcV) - 10 10 for 50% 
mortality of human beings. The situation could be worse because of various leptonic and hadronic radiation 
pathways in the GRB leading to significant fluxes of > 100 MeV - TeV radiation and enhanced cascade UV 
flux (Dar & De Rujula 2002). 

Deflection of charged particles by the Galactic magnetic field means that cosmic-ray protons arrive 
after the delay time 8t ~ 2r^(9 — sin 8) ~ r^9 3 /3c for 9 <C 1, where the characteristic Larmor radius 
tl = O.^s/B^G kpc is assumed to be much larger than the source distance, and 9 is a characteristic 
deflection angle. For a GRB 1 kpc away, >4x 10 18 cV protons arrive within angle 9 of the source direction 
over « 1OOO0 3 years, taking B^q « 4. Diffusive scattering spreads out the arrival times and directions 
compared to this estimate, and implies a "chirping" behavior where higher energy particles arrive first. 
From the preceding expression, the energy dependence of the time delay for a proton with Lorentz factor 
10 9 79 is given by At £* 10 2 d 3 kpc B^ G / l9 s for 79 > B^ G d kpc , and At S 10 2 d 2 kpc / l9 s for 10~ 2 < 79 < 1 , using 
the form of A( 7 ) from Wick et al. (2004). 

Fig. 3 shows the timc-dcpcndcncc of the flux of cosmic-ray neutrons and protons with Lorentz factors 
10 8 < 7 < 10 9 received from a GRB source located one kpc from the Earth. The received flux of prompt 
neutrons with this range of Lorentz factors would be dN n /dAdt 1250/[(6» i /0.1) 2 t dMr (s)] cm s , where 
tdur{s) is the mean duration of the cosmic-ray emission event. Because of the numerical method, this duration 
has been artificially increased to 6.1 x 10 10 s, so that if the actual duration of the event were 600 s, the 
flux would be 10 8 x greater and the duration 10 8 x shorter than shown in the figure, giving constant total 
fluence. The flux units shown in Fig. 3 apply directly to the cosmic-ray protons and neutron-decay protons, 
and we see that the cosmic-ray fluence in the prompt and extended phases are roughly equal. The measured 
energy density of 10 17 - 10 18 eV cosmic rays is w 10~ 18 ergs cm~ 3 , so Fig. 3 shows that the cosmic ray 
energy density at these energies would be « 4 orders of magnitude larger than the currently measured energy 
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Fig. 3. — Fluxes of cosmic rays with 10 < 7 < 10 received at Earth from an on-axis GRB that occurred 1 
kpc away. The fluxes of cosmic-ray neutrons, neutron-decay protons, and protons are shown by the curves 
that decline, rise, and peak during the initial phase. The prompt flux actually lasts for some minutes to 
hours, rather than the 2000 yrs indicated by the numerical simulation, which reflects the approximations of 
the numerical method. 



density during a period of ~ 10 4 yrs following such an event. 

By depositing larger cosmic ray fluence during a much shorter interval, the prompt flux has, however, 
much greater lethality and effect. The occurrence of the sequence of extinction events in the Ordovician due 
to a long lasting ~ 1 Myr ice age (Melott et al. 2004) could happen if the prompt blast induced a long term 
change in the climate, or if the delayed cosmic rays induce a glaciation (Shaviv 2003). On-axis events are 
considerably more damaging than off-axis events which, though more numerous, release < 10 17 eV cosmic 
rays that slowly diffuse towards Earth over periods of thousands of years and longer. 

The hypothesis (Biermann et al. 2004) that the > 10 18 eV cosmic-ray excesses detected with the AGASA 
and SUGAR arrays (Hayashida et al. 1999) are cosmic-ray neutrons from a GRB is not, however, supported 
by our simulations. The relativistic blast waves in a GRB accelerate the highest energy cosmic rays over 
timescales of weeks or less (Zhang & Meszaros 2004), and the high-energy neutrons therefore arrive on this 
same timescale. Cosmic ray protons with energies ~ 10 19 eV are delayed over a timescale « 100006> 3 /£? m g 
years from a source at the distance of the Galactic Center. For the SUGAR excess, which is coincident on 
the sub-degree (9 < 0.02) angular scale with a point source, a GRB would have to take place within weeks of 
the observation for cosmic-ray protons to maintain their direction to the source. Including the requirement 
that the GRB jet was also pointed towards Earth means that an impulsive GRB origin is excluded because 
such an event is highly improbable. The greater (« 10°) extent of the AGASA excess does not conclusively 
exclude a GRB origin, but here the diffuse excess could simply reflect the greater pathlength for cosmic-ray 
proton collisions with spiral arm gas along the Cygnus arm. 

Because the SUGAR point source does not admit an impulsive GRB solution, only cosmic rays from a 
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persistent source, such as a microquasar, could make such an excess. The hypothesis (Dermer 2002; Wick 
ct al. 2004) that GRBs are sources of > 10 14 eV cosmic rays is therefore incompatible with such a source. 
This cosmic-ray origin hypothesis will soon be tested by results from the Auger Observatory 1 to confirm this 
source. If the source is real, then the GRB/cosmic-ray model is incomplete. Searches for neutron /3-decay 
radiation in recent galactic GRBs (Ioka et al. 2004) and around galaxies that host GRBs (Dermer 2002) 
provide further tests of the hypothesis that high-energy cosmic rays are accelerated by GRBs. 

The Earth resides on the inner edge of a spiral arm, and not in an OB association where high-mass 
stars, and therefore GRBs, are usually found. The greater likelihood for intense irradiation events by GRBs 
excludes OB associations from the galactic habitable zone (Lineweaver et al. 2004) except for planets with 
very thick (> 1000 gm cm~ 2 ) atmospheres (Smith, Scalo, & Wheeler 2004) needed to limit the radiation 
effects. 
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